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(54) Method of selecting radio wave propagation paths in a spread spectrum communication 
system 



(57) The signal received in a CDMA system is sam- 
pled at a rate greater than the chip rate of the spreading 
sequences In order to assign the propagation delays to 
the arms of a rake receiver, consecutive delays are test- 
ed within a time window by evaluating, for each delay 
tested, an energy of the correlation between the signal 
received and a reference spreading sequence to which 
said tested delay is applied A selection threshold is de- 
termined from the largest of the energies evaluated. A 
first list is compiled containing delays from the window 
for which the histogram of the evaluated energies ex- 
hibits a local maximum greater than the selection 
threshold, and a second list is compiled containing de- 
lays neighbouring the delays of the first list The delays 
assigned to the arms of the rake receiver are selected 
initially from the first list and then, if this is insufficient, 
from the second list 
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Description 



The present invention relates to the field of direct sequence spread spectrum (DSSS) radio communications. It 
applies in particular to code division multiple access (CDMA) radio telephony systems 

5 CDMA is a mode of digital communication with spread spectrum in which various transmission channels are created 

by using, for each channel, spreading sequences which modulate the information bits to be transmitted The spreading 
sequences have a rate greater than that of the information bits in order to spread the spectrum of the radio signal. 
Their auto- and cross-correlation properties are made suitable to allow the multiplexing of the various channels: they 
are in general mutually orthogonal or quasi-orthogonal, pseudo-random sequences taking the values -1 or +1. These 

io sequences are often complex, that is to say they comprise two real sequences which modulate the same binary stream 
and which are modulated by two quadrature carrier waves. 

The use of CDMA in the field of cellular radio telephony is described in Chapter I of the work "Mobile radio com- 
munications" by Raymond Steele, Pentech Press, London 1 992, and also in the article "On the system design aspects 
of code division multiple access (CDMA) applied to digital cellular and personal communications networks", by A. 

'5 Salmasi and K.S. Gilhousen, Proc. of the 41st IEEE Vehicular Technology Conference, St Louis, MO, 19-22 May 1 991 , 
and in patent US 5 103 459. 

For the downlink, the multiplexed transmission channels are formed at the level of the base station of each cell of 
the network. Each mobile station located in the cell uses a particular spreading sequence in order to retrieve, from the 
overall radio signal emitted by the base station, the information bits which are intended for it. One of the CDMA channels 

20 of the down link is a pilot channel on which the base station emits a reference sequence with no modulation of infor- 
mation bits Reception on the pilot channel makes it possible for the mobile stations to synchronize with one another 
and to evaluate the response of the propagation channel so as to ensure coherent demodulation on the other channels. 

For the uplink, each mobile station emits a radio signal formed using a particular spreading sequence. Correlation 
with this sequence enables the base station to retrieve the information bits transmitted by the mobile station, from the 

25 radio signals received from various mobile stations in the cell. There is no pilot channel in the signal emitted by a mobile 
station. 

For most radio communication systems, the multiple propagation paths between the emitting and receiving stations 
constitute a tricky problem, in particular because they give rise to signal fading (Rayleigh fading). An important advan- 
tage of CDMA is that it lessens the drawbacks related to the multiple paths, by virtue of the spreading of the signal 
30 spectrum. 

With CDMA, it is even possible to profit from the existence of the multiple paths in order to improve reception 
performance, by space diversity techniques. A rake receiver is then used, each arm of which provides for the reception 
of the same signals along a selected propagation path. Each path is identified by a delay which is applied to the 
spreading sequence in order to determine a correlation between the signal received and the delayed sequence. The 

35 correlations obtained in various arms of the receiver can then be combined to restore the information bits transmitted. 

This space diversity technique is explained in the patent US 5109390. In the system described in this document, 
one of the arms of the rake receiver serves to scan the time domain in order to search for new propagation paths. The 
purpose of this search is to select in respect of reception the paths/delays for which the energy of the signal received 
after correlation is greatest. The delays are tested sequentially within a predefined time window by determining their 

40 associated energy. Each tested delay corresponds to a shift by an integer or fractional number of samples or "chips" 
of the spreading sequence which is applied to the latter in order to compute the correlation. The time window of the 
tested delays typically has a duration of the order of 100 jis, this representing a propagation over around 30 km Once 
a delay has been selected and assigned to an arm of the rake receiver, the corresponding path is tracked by means 
of a tracking loop, of the Costas loop type, included within this arm. An example of such a loop is described in the 

->s article "Theory of Spread-Spectrum Communications - A Tutorial" by R.L. Pickholtz et al., IEEE Trans, on Communi- 
cations, Vol. COM-30, No. 5, May 1982. 

In the histogram of the evaluated energies versus the delays tested in the time window, each propagation path 
appears as a pulse whose temporal width is proportional to the inverse of the width of the transmission passband. To 
improve the resolution in the estimation of the delays, the baseband received signal is oversampled, that is to say is 

50 digitized at a rate f s = 1/T B which is a multiple of the chip rate f c = 1/T C of the spreading sequences : f s = M.f c . The main 
lobe of the impulse response of the overall transmission system, assuming a single propagation path, then has a width 
of several samples When a first delay is selected from the main lobe, there is a fairly high probability of also selecting 
one or more additional delays from the M delays preceding and the M delays following this first delay in the time window. 
These additional delays afford a gain in space diversity only if the multiple propagation paths have a concentrated 

55 distribution. Otherwise, these additional delays certainly contribute to increasing the energy received, but they afford 
no gain in diversity. 

An object of the present invention is to propose a method of selecting the delays for a rake receiver which makes 
it possible efficiently to separate the contributions of the propagation paths and those of the noise so as fully to exploit 
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the space/time diversity without increasing the noise in the demodulated output 

The invention thus proposes a method of selecting multiple propagation delays for the reception of messages 
transmitted in a spread spectrum radio communication system, each selected delay being applied to at least one 
spreading sequence in order to determine a correlation between a received signal and said spreading sequence in a 
5 reception arm of a K-arm rake receiver. The received signal is sampled at a sampling rate greater than the chip rate 
of the spreading sequences. Consecutive delays which are multiple of the inverse of the sampling rate are tested within 
a time window, by evaluating, for each tested delay, an energy o( the correlation between the received signal and a 
reference spreading sequence to which said tested delay is applied, and the delays are selected on the basis of the 
distribution of the energies evaluated in respect of the tested delays. The method comprises steps la/ to /d/ below: 

to 

/a/ determining a selection threshold from the largest of the energies evaluated in respect of the delays tested 
within the time window; 

/b/ constructing a first list of delays and a second list of delays, the first list containing delays from the window for 
which the histogram of the evaluated energies versus the delays exhibits a local maximum greater than the se- 
15 lection threshold, and the second list containing delays which are each neighbours of at least one delay from the 

first list: 

Id if the number of delays in the first list is equal to or greater than the number K of arms in the rake receiver, 
selecting the K delays from the first list for which the evaluated energies are largest: 

Idl if the number of delays in the first list is less than the number of arms in the rake receiver selecting the delays 
20 from the first list and at least one delay from the second list for which the evaluated energy is greater than the 

selection threshold. 



In this way, the selection procedure favours the delays in the first list, which correspond in general to distinct 
propagation paths. The second list is ultimately explored, in order to increase the received power without necessarily 
25 contributing to reception diversity, only if the number of delays selected from the first list is insufficient for the K arms 
of the rake receiver. 

The diversity of the multiple paths can thus be better exploited, particularly when the number of arms in the rake 
receiver is limited by reason of restrictions on the hardware complexity and/or energy consumption of the receiver. 
Other features and advantages of the invention will emerge on reading the description below of a non-limiting 
30 embodiment, with reference to the appended drawings in which: 

Figure 1 is a diagram of a CDMA receiver implementing the invention: 

Figure 2 is a diagram of an estimation unit of the propagation channel of the receiver of Figure 1 : 
Figures 3 to 10 are histograms illustrating a method of selecting delays according to the invention: and 
35 - Figure 11 is a flowchart illustrating a way of executing such method. 

In a cellular radio telephony system, base stations are distributed over the territory to be covered In the region, 
or cell, covered by a base station, several mobile stations may communicate simultaneously with the base station. The 
process according to the invention is described below in its application to the downlinks, that is to say to the transmission 
->o of the signals from a base station to mobile stations, for which code division multiple access (CDMA) is used. 

Each base station sets up CDMA transmission channels each defined by a respective spreading sequence CP, 
C1,. .., Cx. The spreading sequences are for example complex, with real and imaginary parts each taking the values 
+1 or -1. These sequences are periodic, pseudo-random and almost uncorrelated. They have a chip rate f c = 1/T C 
greater than the bit rate f b = 1/T b of the information bits to be transmitted, for example a rate of 5 MHz. In each CDMA 
-*5 voice communication channel, the information bits to be transmitted, which represent voice or data signals previously 
coded by conventional techniques, are multiplied by the spreading sequence C1, Cx produced by an appropriate 
generator. Each of the sequences C1 , . Cx is characteristic of the user (mobile station) with which the communication 
is performed over the corresponding channel. 

One of the CDMA channels formed by the base station is a pilot channel, on which a reference spreading sequence 
50 CP is emitted, with no modulation of information bits. 

The baseband spread spectrum signals formed in the various CDMA channels are combined and then modulated 
on a carrier frequency. Combination consists in a summation, possibly weighted, of the baseband signals. The mod- 
ulation can be bi-phase or quadrature phase shift keying (BPSK or QPSK). The carrier frequency is for example 2.4 
GHz. The radio signal obtained after modulation is emitted by the base station towards the mobile stations located in 
55 the cell. 

A receiver of a mobile station is represented in Figure 1 The radio signal received on the antenna 6 is firstly 
processed by the radio stage 7 to yield a baseband complex signal. This complex signal is subjected to a low-pass 
filtering 8 and then digitized by the analogue/digital converter 9 at a sampling rate f s greater than the chip rate f c of the 
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spreading sequences: f 0 = M.f c , with M > 1 , M generally being an integer, typically M = 2 or 4. The CDMA components 
emitted by the base station on various CDMA channels are included in the output signal r(t) from the converter 9 The 
mobile station can retrieve the information bits intended for it by calculating the correlations between the signal r and 
the spreading sequence Cx specific to it. 

5 The receiver of the mobile station comprises a rake receiver 10. In the example represented, the rake receiver 10 

contains K = 2 data reception arms 1 , 2 Each reception arm 1 , 2 provides (or the reception of the signal along a 
propagation path identified by a delay d(1), d(2). The delays d(1 ), d(2) are assigned to the reception arms 1 , 2 by an 
estimation unit 40 of the propagation channel. The rake receiver 10 moreover includes a pseudo-random generator 4 
which delivers the receiver-specific spreading sequence Cx to the arms 1, 2, and a combining unit 5 which receives 

'0 the signals B1 , 82 delivered by the reception arms so as to deduce therefrom a weighted estimation B of the information 
bits received. 

A time base 30 supplies the elements of the receiver with various clock signals obtained by frequency division 
using a crystal oscillator, namely: 

'5 - CK S , sampling clock of period T s supplied to the analogue/digital converter 9: 

CK C , chip clock giving the chip rate of the spreading sequences, of period T c = M.T S ; 
CK b , bit clock giving the bit rate of the information bits, of period T b = M b .T s : 

CK|, clock supplied to the propagation channel estimation unit 40, of period T, = Mj.T s , fixing the integration time 
Tj used to estimate the complex impulse response of the propagation channel, and the updating period (or the 
complex amplitudes A(k)ei'M k) supplied to the arms of the rake receiver in order to provide for coherent demodu- 
lation: 

CK a , clock supplied to the unit 40, of period T a = M a .T s , fixing the update period for the propagation delays d(k) 
supplied to the arms of the rake receiver. 

2S Each arm 1 , 2 of the rake receiver comprises a correlator formed, in the example represented, by a buffer memory 

11, 21 , a complex multiplier 1 2, 22 and a summator-accumulator 1 3, 23. The buffer memory 11, 21 is sized so as to 
receive P successive samples of the received signal r(t), P.T a being the range of possible variation of the propagation 
delays, which is defined as a function of the maximum size of the cells of the network. For example, for a sampling 
frequency f s of 1 0 MHz, we may typically take P - 1024, namely a maximum delay of the order of 100 (.is corresponding 

30 to a propagation path of the order of 30 km. In write mode, the memory 11 , 21 behaves like a shift register. A read is 
performed at each cycle of the clock CK S , at an address which depends on the delay value d(k) (expressed in units of 
the period T s ). In each sampling period, the memory 11, 21 produces the value of the received signal delayed by P-d 
(k) sample durations T s . This value is multiplied by the sample, delivered by the generator 4, of the spreading sequence 
Cx. The initial position of the generator 4 is set as a function of synchronization data received from the base station 

35 on a synchronization channel which is dual of the pilot channel (see for example EP-A-0 630 120), and by taking into 
account a constant delay of P/M chips, so that the delay between the spreading sequence and the signal with which 
it is correlated in the arm k is properly (P/M)T C - [P-d(k)]T s = d(k).T s . 

The output from the multiplier 12, 22 is integrated by the summator-accumulator 13, 23 at the sampling rate f 0 . 
The integration time corresponds to the duration T b of an information bit, supplied by the bit clock CK b A multiplier 14,- 

-*0 24 multiplies the output from the summator 13, 23 by the complex conjugate, supplied by the unit 40, of the complex 
amplitude A(k)ei'i'< k ) of the response of the propagation channel along the path corresponding to the delay d(k). At each 
bit time, the output from the multiplier 14, 24 represents an estimation B1 . B2 of the bit value received along the relevant 
path. These estimations are combined by the unit 5 in order to build the overall estimation B. 

Figure 2 shows an estimation unit 40 of the propagation channel, usable in the receiver of Figure 1. This unit 40 

-*5 comprises a sliding correlator 50 which determines the correlations between the signal received r(t) and the reference 
sequence CP affected by a delay varying in an N-delay search time window enclosed within the P-delay range. In the 
example represented, the sliding correlator 50 includes a buffer memory 51 iden-tical to the memories 11, 21 of the 
arms of the rake receiver 1 0, and a matched filter 53 with a length of M| samples. A counter 43 over log 2 N bits (here 
it is assumed that N is a power of 2, for example N = 1 28) delivers the index n of the tested delay, incremented by one 

so unit in each period of the clock CK S (0<n<N) This index n is added, in 44, to an integer n Q determining the position of 
the smallest delay n Q .T s in the search window, in order to obtain the sum d = n D + n representing the tested delay. The 
integer n 0 is supplied to the memory 51 so that the latter may deliver the sample of the received signal r(t) with a delay 
of P-n D samples. 

The reference sequence CP is supplied to the matched filter 53 by the generator 45 at the chip rate given by the 
55 clock CK C , with a fixed delay of P/M chips. With each overflow of the counter 43, the latter delivers a pulse CfCVN to 
the matched filter 53 in order to update the coefficients of the matched filter 53 Thus, during the ith cycle of the signal 
CK S /N of period NT S , that is to say between the instants t=iNT s and t=(i+1 )NT S . the M, coefficients of the matched filter 
53 are the samples of the reference sequence CP at the sampling instants (iN-P-Mj+1 )T S to (iN-P)T s . At the instant t= 
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(iN+n)T s , the matched filter delivers the correlation: 



r[t-{P-a*j) T ) ] *CP[ (lN-P-j) Tj 

O B 0 



= £ r(t'-jT g ) *CP(t'-jT g -dT g ) 
1=0 



with t' = T-(P-n Q )T 3 . The delay d=n 0 +n is therefore properly tested when the counter 43 displays the value n, by calcu- 
lating a correlation between the input signal r and the pilot sequence CP over an integration time Tj = Mj.T s . 

The correlation delivered by the sliding correlator 50 when the counter 43 displays the value n is an estimation of 

20 the complex amplitude A(n)e'* | '' n) of the response of the propagation channel for the propagation delay of (n 0 + n)T s . 
This complex amplitude is stored at the address d = n Q + n of a memory included within a unit 46 for processing the 
complex impulse response (CIR). possibly after suitable processing (filtering, conjugation, multiplication by a scale 
factor). The unit 46 knows the relative delays d(k) applied to the arms of the rake receiver 10, and is thus able to read 
from its memory at the addresses d(k) the values of the complex amplitudes supplied to the multipliers 14, 24 of the 

25 arms of the receiver. The frequency of updating the complex amplitudes is, in the case illustrated in Figure 2 in which 
the delays are tested by means of a matched filter, equal to the frequency of the clock CK, giving the integration time Tj. 

The output from the multiplier 53 is moreover supplied to an element 47 which calculates its modulus squared and 
supplies it to a unit 48 for calculating mean values In each period of the clock CK S , the unit 48 receives the relative 
value n of the tested delay and an estimation w(n) of the energy of the correlation between the signal received and the 

30 reference sequence affected by this delay. This estimation w(n) is the modulus squared supplied by the element 47, 
the duration of estimation of the energy being the period Tj of the clock CKj. For each value of n (0 < n < N), the unit 
48 accumulates several successive values of the estimation w(n) in order to obtain an average W(n) thereof. The 
accumulation variables are reset to zero at the beginning of each cycle of the clock CK a . At the end of each cycle of 
the clock CK a , the unit 48 for calculating mean values makes the averages W(n) available to a path selection unit 49. 

•35 The average W(n) is an evaluation of the energy received in respect of the relative delay n, obtained by averaging M a / 
Mj successive values of the estimation W(n). It is on the basis of the evaluated energies W(n) that the unit 49 selects 
the delays d(k). Once these delays have been selected, the unit 49 supplies them to the arms of the rake receiver 10 
and to the unit 46 for processing the CIR 

The propagation delays d(k) are in general stable over longer times than the complex amplitudes A(k)e)'t' (k >. The 

■to estimation unit 40 of Figure 2 allows this to be taken into account by retaining an update period J, for the complex 
amplitudes which is shorter than the update period T a for the delays, for example T a ~ 10 ms and Tj — 250 ms. corre- 
sponding to M a /Mj ~ 40 estimations w(n) taken into account for each evaluated energy W(n) 

Figures 3 to 9 illustrate a procedure for selecting the delays on the basis of the energies W(n) evaluated during 
an update period of duration T a in a particular case where K = 4. M = 2, N = 128. Figure 3 shows an example of a 

-> s histogram of W(n) supplied to the path selection unit 49 by the unit 48 for calculating averages. A selection threshold 
TH is firstly defined by deducting a margin A from the largest value of W(n), expressed in decibels. The delay n = 25 
for which the energy W(n) is largest is included in a first initially empty list L1 . M = 2 neighbouring delays (n = 24 and 
n = 26) are included in a second initially empty list L2. 

A propagation path is manifested by a peak X1 -X5 in the histogram of the energies, the width of which corresponds 

so toM + 1 samples, or to M samples in the particular case in which the delay is very close to an odd integer multiple of 
half the sampling period T s (the case of the peaks X3 and X4 in Figure 3). In the particular case where M = 2, a peak 
of width of M + 1 =3 generally exhibits a maximum for the central sample n (coincidence between the two chip values 
over the whole duration of the sample) and a relatively high value again for the samples n - 1 and n + 1 (coincidence 
between the two chip values over part only of the duration of the sample). It is therefore appropriate to take the delays 

55 n - 1 and n + 1 as the two "neighbours" to be included in the list L2. If the peak has a width M =2 only, one of the two 
delays included in the list L2 represents a low-energy correlation, but this delay will ultimately have a low probability 
of being selected. 

If M > 2, the chip values coincide over M - 1 > 2 consecutive sample durations, so that the energy maximum does 
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not necessarily correspond to the centre of the peak. One possibility is to retain as "neighbours" the M/2 delays pre- 
ceding and the M/2 delays following the delay for which the energy is a maximum. It is also possible to analyse the 
shape of the peak so that more neighbours may be taken from one side of the maximum than from the other in order 
to take account of the cases in which the maximum in the peak is not central. In all cases, a delay "neighbouring" 
£ another delay is at most M-1 samples distant from this other delay. 

In an alternative embodiment, a neighbouring delay is included in the second list L2 only on condition that the 
energy evaluated in respect of this delay is greater than the threshold TH. This does not change the delays which will 
ultimately be selected since the energy evaluated in respect of a selected delay should always be greater than the 
threshold. 

io Once the first delay has been included in the list L1 and its neighbours in the list L2, these M + 1 delays are deleted 

from the histogram n,W(n) (Figure 4). We then search for the remaining delay for which the energy W(n) is a maximum 
(n = 47 in the example represented). If this maximum is greater than the threshold TH, this delay is included in the list 
LI and its M = 2 neighbours (n = 46 and n = 48) are included in the list L2: 

The lists construction procedure mentioned in the previous paragraph is repeated until the histogram n, W(n) no 

15 longer exhibits any energy greater than the threshold TH (Figure 7). Figures 5, 6 and 7 show the histogram n, W(n) 
after 2, 4 and 5 iterations respectively. Figure 8 shows the ultimate contents of the second list L2, the 10 delays being 
arranged in ascending order. 

The selection of the delays is firstly performed from the first list L1 . If this list L1 does not contain sufficient delays 
for the K arms of the rake receiver, the second list L2 is examined so as to select other delays therefrom. 

20 in the case of Figure 9, the rake receiver includes K - 4 arms With the list L1 containing 5 delays, the first 4, n(l ) 

to n(4), are selected, that is to say those for which the energies W(n) are largest 

In the case of Figure 10, the rake receiver includes K = 8 arms. The 5 delays n(1) to n(5) from the list L1 are 
selected, as are the 3 delays n(6) to n(8) from the list L2 for which the energies W(n) are largest 

The absolute delays d(k) ultimately assigned by the unit 49 to the arms of the receiver are the relative delays n(k) 

2S increased by the position n 0 of the first sample from the search window. 

Note that, in the case of Figure 10, the selected delays n(1 ) to n(8) correspond to the 8 delays for which the energies 
are the largest at the start (Figure 3). In this particular case the selection procedure leads to the same result as a 
conventional procedure based only on the energy levels. However, in the case of Figure 9, where the number of arms 
is smaller, 4 delays are selected corresponding to 4 distinct paths, whereas the conventional procedure would lead to 

oo the selection of the delays n = 24, 25, 26 and 47, including 3 neighbouring delays which are probably due to the same 
propagation path. The invention therefore makes it possible to derive greater benefit from the diversity in reception 
when the number of arms is limited. 

Figure 11 shows an example of a flowchart for the procedure for selecting the delays from the histogram n,W(n) 
obtained over a period T a . In the first step 60, the unit 49 calculates the threshold TH by deducting the margin A from 

os the largest of the energies W(n) over the window, and initializes the integer i to 0. A loop 62. 64, 66 is then executed 
several times to construct the initially empty lists L1 and L2. In step 62, the delay j is determined for which the evaluated 
energy is the largest in the search window with the exclusion of the delays in the lists L1 and L2. This largest energy 
is compared with the threshold TH in step 64. If the comparison 64 shows that W(j) > TH, step 66 is executed before 
returning to step 62. In step 66, the integer i is incremented by one unit and the delay n(i) is taken equal to the integer 

-*o j which previously maximized the energy. The delay n(i) is included in the list LI, and the neighbouring delays (n(i)-l 
and n(i)+1 in the case where M = 2) are included in the list L2. 

When the comparison 64 shows that W(j) < TH, the loop 62, 64, 66 is exited, the lists L1 and L2 then being complete. 
The number i of delays in the first list is compared, at 70, with the number K of arms in the rake receiver. If i > K, the 
first K delays in the first list L1 are selected in step 72, the absolute delays d(k) being obtained for 1 < k < K from d(k) 

■*5 = n Q + n(k). If the first comparison 70 shows that i < K, the first list L1 does not contain sufficient delays for the K arms 
of the receiver, and the second list L2 is processed in a loop 74, 76, 78, 70. In step 74 the delay j is determined for 
which the evaluated energy is the largest in the list L2, with the exclusion of an initially empty list L3. In step 76, the 
maximized energy W(j) is compared with the threshold TH If W(j) > TH, step 78 is executed before returning to the 
comparison step 70. In step 78, the integer i is incremented by one unit, and the delay n(i) is taken equal to the integer 

so j which previously maximized the energy. The delay n(i) is included in the list L3. When the comparison 70 shows that 
i = K, K delays are selected in step 72, namely all the delays in the list L1 and one or more delays in the list L2 (those 
for which the energies are the largest). If an instance of- the comparison 76 shows that W(j) < TH, the second list L2 
does not contain sufficient delays greater than the threshold TH to supply a total of K delays together with those from 
the first list. In this case, the delays are selected in step 80, equal in number to the current value of the integer i. 

55 Once the delays have been selected in step 72 or 80, the unit 49 calculates two quantities E s and l Q in step 82 

and then determines, in step 84, the parameters A and n c for the following period T a . The new position n 0 of the search 
window is determined from the distribution of the delays which have just been selected. A simple possibility is to take 
n Q equal to the smallest of the delays d(k), less a constant value. More sophisticated algorithms can, however, be 
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applied in order to position the window. 

The margin A is determined on the basis of estimations 'of the signal-to-interferers ratio (SIR). The overall energy 
received E s has been determined in step 82 by summing the energies W[n(k)j evaluated in respect of the delays n(k) 
retained in step 72 or 80, whereas the power density l 0 of the interf erers has been determined by summing the energies 
s W(n) evaluated in respect of those delays in the search window which are not included in the lists L1 and L2. The 
quotient EJ\ 0 represents an estimation of the signal-to-interferers ratio in each period T a . These estimations are sub- 
jected to low-pass filtering, and the margin A is taken equal to the filtered estimation E s /I c multiplied by a predetermined 
coefficient. 

By way of illustration, the applicant has tested the process in a particular case in which the sampling rate f s was 
'0 10 MHz, for a chip rate of 5 Mchips/s (M = 2). The search window covered N = 1 28 samples, i.e. 1 2.8 us. In the sliding 
correlator, the integration time T; corresponded to Mj = 2560 samples or 1 280 chips. The update period for the complex 
amplitudes was equal to the integration time T h i.e. 256 (.is. The update period T a for the delays was 10.24 ms (M a = 
102400). M a /Mj = 40 estimations w(n) were taken into account in evaluating the energies W(n). With a receiver with K 
~ 4 arms, the applicant has observed better exploitation of the delays corresponding to paths other than that having 
'5 the highest energy, in contrast to the conventional receivers which have too great a tendency simultaneously to select 
several neighbouring delays corresponding to the path of highest energy. 

The invention has been described above in an illustrative application to the downlinks of a radio telephone system. 
It can of course be applied in similar fashion to the uplinks when the receiver of the base station includes a diversity 
receiver (see US 5109390). In this case, the reference sequence used in the calculation of the reception energies is 
20 no longer a sequence emitted over a pilot channel, but rather the spreading sequence serving for the communications 



Claims 

25 1 . A method of selecting multiple propagation delays (or the reception ol messages transmitted in a spread spectrum 
radio communication system, each selected delay (d(1),d(2)) being applied to at least one spreading sequence 
(Cx) in order to determine a correlation between a received signal and said spreading sequence in a reception 
arm (1,2) of a K-arm rake receiver (10), in which method the received signal is sampled at a sampling rate (Ts) 
greater than the chip rate (Tc) of the spreading sequences, consecutive delays which are multiple of the inverse 

30 of the sampling rate are tested within a time window, by evaluating, for each tested delay, an energy (W(n)) of the 

correlation between the received signal and a reference spreading sequence (CP) to which said tested delay is 
applied, and the delays are selected on the basis of the distribution of the energies evaluated in respect of the 
tested delays, characterized in that the selection of the delays comprises the steps of: 

35 /a/ determining a selection threshold (TH) from the largest of the energies evaluated in respect of the delays 

tested within the time window: 

/b/ constructing a first list of delays (L1 ) and a second list of delays (L2), the first list containing delays from 
the window for which the histogram of the evaluated energies versus the delays exhibits a local maximum 
greater than the selection threshold, and the second list containing delays which are each neighbours of at 

->o least one delay from the first list: 

Id if the number of delays in the first list (L1) is equal to or greater than the number K of arms in the rake 
receiver (10), selecting the K delays from the first list for which the evaluated energies are largest: 
/d/ if the number of delays in the first list (L1 ) is less than the number of arms in the rake receiver (10), selecting 
the delays from the first list and at least one delay from the second list (L2) for which the evaluated energy is 

■*5 greater than the selection threshold (TH). 

2. A method according to Claim 1 characterized in that the sampling rate (Ts) is equal to M times the chip rate (Tc) 
of the spreading sequences, M being an integer larger than 1, and in that each delay in the second list (L2) is at 
most M-1 samples distant from a delay in the first list (L1 ). 

so 

3. A method according to Claim 2, characterized in that step /b/ of constructing the lists comprises, starting from 
initially empty lists, several iterations of the succession of the following sub-steps: 

/bl/ determining the tested delay, from among the delays in the time window with the exclusion of the delays 
55 in the first and second lists (Li, L2), for which the evaluated energy (W(n)) is largest: 

/b2/ comparing the evaluated energy in respect of the delay determined in sub-step /bl/ with the selection 
threshold (TH): and 

/b3/ if sub-step /b2/ shows an energy greater than the selection threshold, including in the first list (L1) the 
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delay determined in sub-step /bl/, and including in the second list (L2) at least each delay, from among M 
delays neighbouring the delay determined in sub-step /bl/, for which the evaluated energy is greater than the 
selection threshold. 

A method according to Claim 3, characterized in thai M being an even integer, said M delays neighbouring a delay 
determined in sub-step /bl/, from among which at least one delay is included in the second list (L2) in the next 
sub-step /b3/. are the M/2 delays preceding and the M/2 delays following said delay determined in sub-step /b1/. 

A method according to any one of Claims 1 to 4, characterized in that the lists (L1 , L2) are constructed in step /bl 
in such a way that none of the energies evaluated in respect of the delays in the time window with the exclusion 
of the delays from the first and second lists is greater than the selection threshold (TH), the method furthermore 
comprising a step of estimating a signal-to-interferers ratio based on a quotient between on the one hand the sum 
of the energies evaluated in respect of the delays selected in step Id or step /d/ and on the other hand the sum 
of the energies evaluated in respect of the delays in the time window with the exclusion of the delays from the first 
and second lists (L1, L2). 

A method according to Claim 5, characterized in that the delays are selected per successive update periods of M a 
samples, and in that the selection threshold (TH) determined in step /a/ for an update period is calculated as a 
function of the largest of the energies evaluated in respect of the delays tested within the time window during said 
update period and of the previously estimated signal-to-interferers ratio. 

A method according to any one of Claims 1 to 6, characterized in that each arm (1,2) of the rake receiver (10) is 
arranged to perform a coherent demodulation in respect of the propagation delay (d(k)) assigned thereto, on the 
basis of an estimation of a complex reception amplitude (A(k)ei'l' (k )) associated with said delay, and in that the 
estimations of the complex amplitudes (A(k)e'*( k) ) are updated at a frequency higher than the frequency of updating 
the delays. 
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